Momentum distribution of a trapped Fermi gas with large scattering length 



L. Viverit, 1,2 S. Giorgini, 1 L. P. Pitaevskii, 1,3 and S. Stringari 1 

1 Dipartimento di Fisica, Universita di Trento and BEC-INFM, 1-38050 Povo, Italy 
2 Dipartimento di Fisica Universita di Milano, via Celoria 16, 1-20122 Milan, Italy 
s ' Kavitza Institute for Physical Problems, 117334 Moscow, Russia 

Using a scattering length parametrization of the BCS-BEC crossover as well as the local density 
approximation for the density profile, we calculate the momentum distribution of a harmonically 
trapped atomic Fermi gas at zero temperature. Various interaction regimes are considered, including 
the BCS phase, the unitarity limit and the molecular regime. We show that the relevant parameter 
which characterizes the crossover is given by the dimensionless combination N 1 ^ 6 a/ai lo , where N is 
the number of atoms, a is the scattering length and au is the oscillator length. The width of the 
momentum distribution is shown to depend in a crucial way on the value and sign of this parameter. 
Our predictions can be relevant for experiments on ultracold atomic Fermi gases near a Feshbach 
' resonance. 

o 

Recent experiments on ultracold atomic Fermi gases near a Feshbach resonance, have pointed out the crucial role 
played by two-body interactions. Spectacular effects concern the hydrodynamic behaviour exhibited by the expansion 
after release of the trap 0, 0, Q and the evidence of molecular formation for positive values of the scattering length 
These experiments open important perspectives towards the realization of the superfluid phase in Fermi gases, 
whose critical temperature has been predicted to be significantly enhanced by resonance effects 

Bui, 

and towards 

the study of the BCS-BEC crossover [l(J, HH H2 and of many-body effects in the presence of a Feshbach resonance 

^: mum 

qq . The purpose of this paper is to point out that an important diagnostics of the state of the system across the 
' resonance is the atomic momentum distribution whose shape, at low temperature, is very sensitive to the size and 
sign of the scattering length. The momentum distribution can be directly measured in these systems either by suddenly 
■ switching off the scattering length and imaging the expanding atomic cloud , or by Bragg spectroscopy 0] . 

In this paper we will consider a Fermi gas of two spin species. The interspecies interaction is characterized by 
the s-wave scattering length a which will be assumed to be larger than the effective range ro of the interaction. 
Furthermore we will also assume that the relative momentum hk between the colliding atoms is smaller than /ro- 
If the scattering length a is large and positive a weakly bound molecule forms in the vacuum with binding energy 
Ei, = —h 2 /ma 2 and size a [17\. A peculiar feature of the resonance is that when k\a\ 3> 1 the scattering amplitude / 
is unitarity limited and / — > i/k, independent of the value of a |18| . 

In the many-body problem the above lengths should be compared with the average distance between particles 
£h ! proportional, at T = 0, to the inverse of the Fermi wavenumber 

kF = (6. 2 ) i/3 «y 3 , (i) 



where the densities n a of the two species will be taken equal. By assuming that both the average distance between 
particles and the modulus of the scattering length are larger than the effective range ro , the effects of interactions will 
not depend on the actual value of ro, but only on the combination kpa. 

For small and negative values of kpa, a homogeneous gas at zero temperature exhibits a BCS superfluid phase, 
whose momentum distribution differs very little from the step function 9(kp — k) of the non- interacting Fermi gas. 

When either the value of the scattering length or of the density increases, correlations become more and more 
important. Eventually, in the unitarity limit kp\a\ 3> 1, the configuration of the system is expected to exhibit a 
universal behavior, independent of the actual value of a |2J|. When we move to the side of the Feshbach resonance 
where the scattering length is large and positive, bound molecules of size a would form in the vacuum. However if 
a is still much larger than the interparticle spacing the actual state of the many-body system will be very different. 
Only when kpa becomes smaller than unity will these molecules not be perturbed by the medium and behave just 
as independent bosons. In this limit the atomic momentum distribution is expected to become ~ l/(fc 2 a 2 + l) 2 , 
corresponding to the momentum distribution of atoms inside molecules. 

In the following we shall work at T — and take into account the effect of the harmonic trapping. This is important 
because the actual values of the densities, and hence of fcp, realized in experiments depend in a crucial way on the 
trapping parameters. We use the model developed in [lj to investigate the BCS-BEC crossover at T = 0, which 
accounts for all the regimes introduced above (BCS limit, unitarity limit, BEC molecular limit) with a unique tuning 
parameter fixed, in uniform gases, by the combination kpa. This approach consists of a generalization of the gap 
equation and the number equation of the usual BCS theory to the whole resonance region, and corresponds to using 
the saddle point approximation for the zero temperature partition function. 
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Eliminating as usual the bare interatomic potential in favor of the scattering t-matrix in the vacuum, one obtains 
the result 
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for the gap equation,where e k = h 2 k 2 /2m, E k = \J (ifc — m) 2 + A 2 \ f 1 1S tne chemical potential, A is the gap parameter 
and V is the volume of the system. The number equation, on the other hand, takes the form 

The two equations form a closed set of equations in the unknowns A and /x and can be solved for a given value of 
kpa. The approach is reliable when kpa — > 0^ with \a\ 3> ro, and approximate in the region kp\a\ > 1. 
It is immediate to see that the particle distribution of atoms of either species is given by the expression 

and corresponds, in the usual Bogoliubov theory of Fermi supcrfluids, to n k = \v k \ 2 where v k are the Bogoliubov 
amplitudes. The above expression for n k is only valid for values of k smaller than the inverse of the effective range 
r . The values of k can however be much larger than 1 /a since we assume a»rj. 

Solutions of equations an d © are available for the whole range of values of kpa [2ll 1221] . Particularly simple 
solutions exist in the limits kpa — > O*. In the BCS limit (kpa — > - ) the solutions are well known: /i ~ ep and 
A ~ &e~ 2 ep exp(7r/2fc^a), with ep = h 2 k 2 F /2m being the Fermi energy. Notice that, since the present theory neglects 
the fluctuations in the number-density and spin-density, the Gorkov correction to the gap A is not included |l9i l2f)j | . 
However, this has no major relevance for the calculation of the particle distribution which in this limit is close to the 
free-gas value. 

In the BEC limit instead one finds /i ~ —h 2 /2ma 2 + 2ttH 2 an a / 'm and A ~ (16/37r) 1 / 2 ei? / ' \Jkpa. The first term in 
the chemical potential is simply Eb/2, i.e. the molecular binding energy per atom. By rewriting the second term as 
2 r Kh 2 aB , n a /rag, where ag = 2a and mg = 2m, one can see that this accounts for the mean field molecule-molecule 
interaction energy (28j . The quantity A instead plays the role of an order parameter for the molecular Bose-Einstein 
condensate (see for example |23|)- 

In both the BCS and BEC limits A is much smaller than the abolute value of the chemical potential, so that the 
particle distribution can be expanded in terms of A/|/i|. In the BCS limit the lowest order in the expansion is simply 
the free Fermi step function, and it is enough to stop at this level for the considerations which follow. In the BEC 
limit one instead finds the result 

nu = ^(kpaf {k2 J +i)2 . (5) 

This is proportional to the square of the Fourier transform of the molecular wavefunction ip(r) = (V 7 *) ex P( — r / a )i 
indicating that in this limit the particle distribution arises from the motion of atoms inside the molecules, and the 
prefactor accounts for the fact that there are n a molecules per unit volume. 

Let us now consider the most interesting unitarity limit where kpa — > ±oo. The left hand side of Eq. @ vanishes, 
and this reduces to an implicit equation for the ratio A//i. A straightforward calculation [2ll |22| gives A//j ~ 1.16 
and from Eq. @ one finally finds the results A ~ 0.69 ep and /i ~ 0.59 ep. In the unitarity limit the proportionality 
of the chemical potential and of A with ep can be deduced from general dimensional arguments and follows from 
the fact that a is no longer a relevant scale in the problem. In this limit the chemical potential is often written in 
the form fj, = (1 + /3) ep where, in our case, (3 ~ —0.41. It is interesting to compare this prediction with the results 
of more elaborated approaches. A recent microscopic calculation based on quantum Monte-Carlo techniques gives 
the value f3 = —0.56 j25j, while the theory of resonance superfluidity based on saddle-point approximation yields the 
value (3 ~ —0.35 j2(J. The value of (3 has also been the object of recent experimental investigations in trapped Fermi 
gases at finite temperature, through the study of the release energy 

The predictions of the theory for the particle distribution of a uniform gas are reported in Fig. ^ which shows 
that during the transition from the ideal gas (dashed line) to the molecular limit (full line) across the unitarity limit 
(long-dashed line) the atomic momentum distribution varies in a dramatic way, reflecting the critical role played by 
two-body interactions. This behavior of the momentum distribution was fi rst p ointed out in |lOllll| . More recently, it 
has been investigated in and, in the unitarity limit, by the authors of |25| . In particular the results of |2^| rather 
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well agree with the curve of Fig. ^for the unitarity limit. The width of the momentum distribution ranges from 
~ kp in the BCS limit to ~ 1/a in the BEC regime. Notice that for the value kpa = 0.5 reported in the figure, the 
momentum distribution calculated by solving Eqs. J2J and (0), is practically indistinguishible from the "molecular" 
formula JSJ). 

In trapped configurations we shall calculate the momentum distribution by introducing the local semiclassical 
particle distribution 

"k(r) = 7T 1 „ O . ( 6 ) 



2 V E k {r) 

where A(r) and /z(r) are obtained by solving Eq. (J2J) and @ with the local value of the density. The momentum 
distribution is then evaluated by integrating the particle distribution in coordinate space: 

»« = / Mr) • (7) 

The momentum distribution coincides with the particle distribution function only in a uniform system, where one 
has the simple relation n(k) = V n^. A fundamental ingredient for the calculation is the density distribution, whose 
shape can vary very much depending on the regime considered and on the trapping parameters [29| . In the following 
the density profile will be determined using a local density approximation, based on the solution of the equation 

H = ^(n a (r)) + V ext (r) . (8) 

In the deep BCS limit the particle distribution is given by the free distribution. For harmonic trapping one finds 
the well known results n<j(r) = n a (0) [1 — (r/Rj'p) 2 } 3 ^ 2 for the density profile, where Rj-p — (48A ? (T ) 1 / 6 ah is the 
Thomas- Fermi radius, and n CT (0) = (4A^ (T /37r 4 ) 1 / 2 aZ° is the central density. N a is the number of fermions in each 
internal state, and a^o = \J fc/mu; the harmonic oscillator length, with u = (ujxUjyUJz) 1 ^ 3 . The momentum distribution 
takes the following form 

where the Fermi wavenumber satisfies the relation 

kpa=(48)^Ny e — , (9) 
a-ho 

and coincides with the value obtained using Eq. Q with the central density n^O). 

In the unitarity limit the density profile has the same form as in the ideal case, since the chemical potential in the 
uniform system has the same power law dependence on the density (fi cx n 2 / 3 ). The only difference with respect to the 
ideal case is an overall rescaling factor: \i = (l + f3)ep , resulting in a contraction of the Thomas-Fermi radius according 
to the law Rj, F — > R TF (1 + /3) 1 / 4 . Similarly, the central density is rescaled according to n(0) — > n(0)(l + /3)~ 3 / 4 and 
consequently the Fermi wavenumber, calculated in the center of the trap, is fixed by 

^-jeg^JL. do) 

(1 + /5) 1 / 4 Oho 

The momentum distribution (J7J| can then be calculated by integrating © , employing the corresponding values for A 
and n and the rescaled density profile. After introducing the dimensionless variable f = r/(l + flY^R^p one finds 
the result: 

„(k) = ^) 3 (^) 3/4 l\rr> (11) 








x 1 



(k/k F f-{l + pf/ 2 f{f) 



V[(W) 2 - (1 + P) 1/2 f(r)} 2 + (1 + /3)"W)(A/eF)] 2 



where /(f) = 1 - f 2 , (3 = -0.41 and A/e F = 0.69. 



4 



Let us finally discuss the molecular BEC limit. To the extent that the condition kpa <g; 1 is satisfied in the center 
of the trap, one can use the molecular distribution function (JSJ everywhere and the resulting momentum distribution 
will be consequently given by 

"( k ) = ° £° {k2 J + 1)2 • (") 

Also in this case it is important to calculate the quantity kpa in terms of the relevant parameters of the trap, and 
to this purpose one has to determine the density profile. The leading density dependent term 2ith 2 an<j/m in the 
chemical potential arises from the molecule-molecule mean-field interaction. This corresponds to the usual interaction 
term in the Gross-Pitaevskii theory for bosons interacting with scattering length as — 2a. In this limit the local 
density approximation yields the T = Thomas-Fermi profile 

n a (i)= n<r (0)(l--^-), (13) 



where 



( 



2/5 



MO) = T ~ ' ( M ) 

4:7raa ho V 2a ho J 

is the central density and Rtf = aho(^5N a a/2ai lo ) 1 ^ 5 is the Thomas-Fermi radius. The use of the local density 
approximation is justified since the Thomas-Fermi parameter N a a/aho is much larger than unity. In deriving the 
above expressions we have assumed that atoms and molecules are trapped with the same oscillator frequency u>. 
Consequently the boson harmonic oscillator length is given by aho/VZ- By evaluating the Fermi momentum at the 
center of the trap one finally finds 

kpa = 2.2 ( Ny e -^— ) . (15) 

V a ho J 

The quantity kpa should be sufficiently small in order to apply result (|12[) for the momentum distribution. It is 
remarkable that in all of the regimes the product kpa, related to the gas parameter na 3 , is fixed by the combination 
Nl - a/ Oh a . This quantity then permits to characterize the various regimes exhibited by the interacting Fermi gas. 
The universality of this combination is a consequence of harmonic trapping, and of the fact that the equation of state 
fi/ep, evaluated in the uniform phase, has been assumed to be a function of kpa. The above discussion also shows 
that the value of kpa is rather insensitive to the value of N a . 

In Fig. |2]we plot the momentum distribution calculated for three different values of nV & a/ aho- The distributions 
are normalized to unity and momenta are expressed in terms of k/k F where k F is the Fermi wavevector of a non- 
interacting gas. Although these results are based on the approximate theory of 0, El El they provide a first 
useful estimate of the momentum distribution which might stimulate new experimental investigations as well as more 
sophisticated and self-consistent theoretical calculations. 

Let us finally discuss the experimental possibilities to measure the momentum distribution and to probe the behavior 
of the system across the resonance. The technique developed in [f| to quickly turn off the scattering length and to 
immediately release the trap is well suited for this purpose. If the magnetic field is switched off in times shorter 
than the inverse of the binding energy H 2 /ma 2 , the molecules are suddenly dissociated and the atoms will expand 
ballistically. If instead the scattering length is switched off adiabatically the system will form deeper molecular states 
characterized by microscopic sizes. The imaging of the expanding atomic cloud provides a direct measurement of the 
momentum distribution of the initial correlated configuration. An important feature emerging from our results (see 
Figs. n an dGJ) is the presence of large k components in the momentum distribution, with the consequent suppression 
of n(k) at small k. These large k components, whose presence is more and more pronounced as one leaves the ideal 
Fermi gas condition kpa — > 0~, affect sizably the release energy of the atomic cloud fixed by the kinetic energy 30]. 

If one instead releases the trap by keeping the scattering length on, the scenario of the expansion will be completely 
different. In this case one expects that the expansion will be governed by the laws of hydrodynamics. This will be 
the case both at zero temperature, where the system is superfluid, and above T c if kp\a\ 3> 1 due to the resonant 
effect in the collisional cross section. Under these conditions the expansion will be anisotropic [3l| if the confining 
trap is not symmetric and the release energy, as well as the shape of the expanding cloud, will be determined by the 
equation of state of the gas. In the unitarity limit the release energy is given by En/N a = (3/8)(l + /3) 1 / 2 (67V cr ) 1 ' /3 ?ia;, 
while in the molecular regime the expansion will be similar to the one of a dilute Bose gas and one predicts En/N a = 
(l/7)(15N a a/a ho ) 2 / 5 faj. 
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So far the discussion has been restricted to zero temperature. An interesting question concerns the behavior of the 
momentum distribution at finite temperatures. For temperatures just above T c , the momentum distribution can be 
determined through the calculation of the single particle Green's function by summing the particle-particle ladder 
diagrams along the lines of |1C( . In particular, for negative scattering length and kp\a\ <C 1, where ksT c <C &f, 
one finds that the momentum distribution of the free Fermi gas is slightly broadened due to thermal effects. In 
the opposite molecular regime (a > and kpa <C 1) one recovers the momentum distribution JSJ) of atoms inside 
molecules, provided the binding energy of the molecule is much larger than the temperature h 2 /ma 2 3> ksT . In fact 
at the critical temperature for the Bose-Einstein condensation of molecules one finds fcsT/|£?s| = 0.22 (kpa) 2 , and 
the thermal dissociation of molecules is negligible if kpa is sufficiently small. 

In conclusion, in this paper we have reported a first calculation of the momentum distribution of a harmonically 
trapped two-component Fermi gas as a function of the scattering length in the BCS-BEC crossover. We expect that 
the large deviations from the ideal Fermi gas distribution can be addressed experimentally in systems close to a 
Feshbach resonance. 

We thank A. Perali for useful comments. 




FIG. 1: Particle distribution of a uniform gas in the BEC (solid line), unitarity (long dashed line) and BCS regime (short 
dashed line). For the BCS regime we have chosen kpa = —0.5, and for the BEC one kpa — 0.5. In the inset the momentum 
distribution in the BEC regime is compared with the molecular distribution JS) (dotted line). 
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k = k/k£ 



FIG. 2: Momentum distribution of a trapped gas in the BEC (solid line), unitarity (long dashed line) and BCS regime (short 

dashed line). In the BEC regime we have chosen N^ 6 a/a^o = 0.26 corresponding to k F a = 0.75 (and k° F a = 0.5), while for 
the BCS regime we have plotted the free fermion distribution. The momentum distributions are multiplied by k 2 to emphasize 
the large k behavior and are normalized so that J d 3 k n(k) — 1. 
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